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ABSTRACT 

Effects of magnetic coupling (MC) of a rotating black hole (BH) with its surrounding 
accretion disk are discussed in detail in the following aspects: (i) The mapping relation 
between the angular coordinate on the BH horizon and the radial coordinate on the disk 
is modified based on a more reasonable configuration of magnetic field, and a condition for 
coexistence of the Blandford-Znajek (BZ) and the MC process is derived, (ii) The transfer 
direction of energy and angular momentum in MC process is described equivalently by the 
co-rotation radius and by the flow of electromagnetic angular momentum and redshifted 
energy, where the latter is based on an assumption that the theory of BH magnetosphere 
is applicable to both the BZ and MC processes, (iii) The profile of the current on the BH 
horizon and that of the current density flowing from the magnetosphere onto the horizon 
are given in terms of the angular coordinate of the horizon. It is shown that the current on 
the BH horizon varies with the latitude of the horizon and is not continuous at the angular 
boundary between the open and closed magnetic field lines, (iv) The MC effects on disk 
radiation are discussed, and a very steep emissivity is produced by MC process, which is 
consistent with the recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy 
MCG-6-30-15 by a variety of parameters of the BH-disk system. 

Subject headings: Black hole — accretion disk — magnetic fields 



1. Introduction 

Recently magnetic coupling of a rotating black hole (BH) with its surrounding disk has been 
investigated by some authors (Blandford 1999; Li 2000; Li & Paczynski 2000; Li 2002a, 2002b,2002c, 
hereafter Li02a, Li02b and Li02c, respectively; Wang, Xiao & Lei 2002, hereafter WXL; Wang, Lei 
& Ma 2003, hereafter WLM), which can be regarded as one of the variants of the Blandford-Znajek 
(BZ) process proposed two decades ago (Blandford & Znajek 1977). Macdonald and Thorne (1982, 
hereafter MT) reformulated and extended the BZ theory in 3+1 split of Kerr spacetime. In MT the 
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transportation of energy and angular momentum from a rotating BH to a remote astrophysical load is 
not only described by the flow of electromagnetic angular momentum and energy but also by a general 
relativistic version of DC electronic circuit theory based on a stationary, axisymmetric magnetosphere 
anchored in the BH and its surrounding disk. In the BZ process energy and angular momentum are 
extracted from a rotating BH, which is connected with the remote astrophysical load by open magnetic 
field lines. Unfortunately the remote load has been known very few, resulting in some uncertainty in 
the BZ process. 

With the closed magnetic field lines connecting a BH with its surrounding disk the rotating energy 
of the BH provides an energy source for disk radiation, and henceforth this energy mechanism is referred 
to as magnetic coupling (MC) process. The load in MC process is accretion disk, which is much better 
understood than the remote load in the BZ process. Since the magnetic field on the BH horizon is 
brought and held by the surrounding magnetized disk, the magnetic connection between the BH and 
the disk is natural and reasonable. As argued in Li02b, the magnetic connection can produce a very 
steep emissivity compared to the standard accretion, and this result is consistent with the recent XMM- 
Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15. 

Considering the astrophysical importance of MC process, we are intended to improve our previous 
model of MC process given in WXL. Based on a reasonable consideration on the configuration of the 
magnetic field connecting a BH with its surrounding disk we modify the mapping relation between the 
angular coordinate on the BH horizon and the radial coordinate on the disk given in WXL (henceforth 
MRWXL), and derive a condition for the coexistence of the BZ and MC processes (henceforth CE- 
BZMC). The transfer of energy and angular momentum in MC process and the profile of the current 
on the horizon are discussed in detail. In addition, a very steep emissivity arising from MC process is 
produced in our model, being in accord with the recent XMM-Newton observation of the nearby bright 
Seyfert 1 galaxy MCG-6-30-15. 

This paper is organized as follows. In Section 2 MRWXL is modified based on more reasonable 
consideration of the magnetic field in BH magnetosphere, and the angular boundary between the open 
and closed field lines on the horizon is determined naturally by the new mapping relation and the 
assumption of precedence of the magnetic flux penetrating the disk. The condition for CEBZMC is 
derived and described in the 2-dimension parameter space consisting of the BH spin and the power law 
index of the variation of the magnetic field on the disk. In Section 3 the transfer direction of energy and 
angular momentum between the BH and the disk is described equivalently by the co-rotation radius 
based on our recent work (WLM) and by the flow of electromagnetic angular momentum and redshifted 
energy given in MT. In addition, the flux of angular momentum and energy is discussed in detail both 
above and below the equatorial plane of the BH. In Section 4 the profile of electric current flowing on 
the horizon and the profile of current density from the magnetosphere onto the horizon are discussed 
based on an equivalent circuit for a unified model of the BZ and MC processes proposed in WXL. In 
Section 5 the MC effects on the disk radiation are investigated, and a very steep emissivity is worked 
out to fit the recent XMM-Newton observation by using the new mapping relation. Finally, in Section 
6, we summarize our main results. 

Although load disk in MC process is much better understood than the remote load in the BZ process, 
a good picture of the magnetic connection between the BH and the disk has not been obtained. In order 
to facilitate the discussion of the MC effects in an analytic way we make the following assumptions: 

(i) Assumption: The theory of a stationary, axisymmetric magnetosphere anchored in the BH and its 
surrounding disk is applicable not only to the BZ process but also to MC process. The magnetosphere 
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is assumed to be force-free outside the BH and the disk. 

(ii) Assumption: The disk is both stable and perfectly conducting, and the closed magnetic field lines 
are frozen in the disk. The disk is thin and Keplerian, lies in the equatorial plane of the BH with the 
inner boundary being at the marginally stable orbit. 

iii) Assumption: The magnetic field is assumed to be constant on the horizon, and to vary as a power 
law with the radial coordinate of the disc. 

iv) Assumption: The magnetic flux connecting a BH with its surrounding disk takes precedence over 
that connecting the BH with the remote load. 

Assumption (iv) is proposed based on two reasons: (i) The magnetic field on the horizon is brought 
and held by the surrounding magnetized disk, (ii) The disk is much nearer to the BH than the remote 
load. 

Throughout this paper the geometric units G — c — 1 are used. 



2. Mapping relation and a condition for CEBZMC 

Recently a model of BH evolution was proposed in WXL by considering CEBZMC, and the config- 
uration of the poloidal magnetic field is shown in Figure 1, where r in and r out are the radii of inner and 
outer boundary of the MC region, respectively. The angle 9m indicates the angular boundary between 
the open and closed field lines on the horizon, and 9l is the lower boundary angle for the closed field 
lines. 

MRWXL was derived based on the conservation of magnetic flux between the horizon and the disk 
with assumptions (i) - (iii). However there are some inconsistencies with MRWXL as follows: 

(i) It is impossible for one closed field line to leave and return the horizon at the same angle 82 = 7r/2 
on the horizon by the symmetry of the magnetic field above and below the equatorial plane of the BH. 

(ii) The boundary condition used in MRWXL says that the normal component of the magnetic 
field on the horizon is equal to that at the inner edge of the disk, i.e. B± = B z {r ms ) . However, it is not 
consistent with the estimation given by numerical simulation: the strength of the magnetic field at the 
horizon is likely greater than that in the disk (Ghosh & Abramowicz 1997 and the references therein). 

(iii) The boundary angle 9m should not be given at random without any reason. 

In this paper MRWXL is replaced by a new mapping relation by considering a more reasonable 
configuration of the magnetic field. Firstly, the angular 9l is assumed less than 7r/2 to avoid inconsis- 
tency of the closed field line above and below the equatorial plane. Secondly, we replace the previous 
boundary condition by 

2nr H B 1 _ = 2nzu D (r ms ) B z (r ms ) , (1) 

where B± and B z (r ms ) are normal components of magnetic field on the horizon and at the inner edge 
of the disk, respectively. In writing equation (1) we have assumed that the two neighboring loops near 
the inner edge of the disk have the same infinitesimal width, resulting in the ratio of B± to B z (r ms ) 
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varying from 1.8 to 3 for < a* < 1. The quantity vod (r ms ) is the cylindrical radius at the inner edge 
of the last stable orbit radius r ms = Mx 2 ms and reads 

w D (r ms ) = M X 2 ms v / l + X m X + 2x™ 6 s ^, (2) 

where the parameters M and a* = J/M 2 are the BH mass and spin, respectively. Thirdly, as shown 
later, the boundary angle 9m will be determined naturally by the new mapping relation with assumption 
(iv) rather than given at random. 

Considering the flux tube consisting of two adjacent magnetic surfaces "i" and "i + 1" as shown in 
Figure 1, we have A^h = A^d by continuum of magnetic flux, i.e., 

B ± 2n (wp) r=rH d6 = -B z 2n (wp j v^) q _ dr, (3) 

where we have w = (E/p) sin 9. The concerning Kerr metric parameters are given as follows (MT): 

S 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 6, p 2 = r 2 + a 2 cos 2 6, A = r 2 + a 2 - 2Mr, (4) 

So we have 

(^Pl„-„ = (SsmO) r , rB = 2A.fr„sin9, (5) 

Following Blandford (1976) we assume that B z varies as 

B z oc r", (7) 

where ^ = r/r ms is defined as a radial parameter in terms of r ms , and n is the power law index indicating 
the degree of concentration of the magnetic field in the central region of the disk. Incorporating equations 
(1) and (7) we have 

B g = B ± [r H /w D ( rma )]C n - (8) 
Substituting equations (4) — (8) into equation (3) we have 



where 



sin 6d6 = — G (a*; £, n) d£, 



(9) 



- 5 - 



G (a*; £, n) = ^xLV^ + alx m 4 s ^ 2 + ^ Xm ^~ 3 (1Q) 

2>/(l + alx^ A s + 2a^-6) (1 - 2 X ^~ 1 + a^-^" 2 ) 

Integrating equation (9) and setting £ = £ in = r in /r ms = 1 at we derive the new mapping relation 
as follows: 

cosd — cosd L = G (a*; £, n) <i£, (11) 
Considering that the closed field line connects 6*m with £ out = r out /r ms , we have 

G(a*;£,n)c&;. (12) 

However £ out and #m cannot be determined simultaneously by equation (12), even if the values of a*, 
n and are given. To determine the boundary angle 9m, we take assumption (iv), stating that the 
magnetic flux of the closed field lines takes precedence over that due to the open field lines. Thus we 
have the curves of ^ out and 6m versus a* with 6l = 0.457T and different values of n as shown in Figure 2. 

From Figure 2 we obtain the following results: 

(i) As shown in Figure 2a, the parameter 6m remains zero for < a* < 1, provided that the power 
law index n is not so great, such as n = 3.0. In this case the parameter £ out remains finite, varying with 
a* non-monotonically and attaining a maximum as a* approaches unity. 

(ii) As shown in Figure 2b and 2c, the parameter 6m remains zero for the smaller value of a*, and 
becomes positive for the greater value of a*, provided that n is great enough, such as n = 4.0, 4.5. In 
this case the parameter £ out approaches infinite as soon as 6m becomes positive. 

The above results arise from the conservation of the magnetic flux and the precedence for the 
magnetic flux connecting the BH with the disk. In our model the magnetic flux sending out from 
the BH horizon should be the sum of and ^l, which are the fluxes connecting the disk and the 
remote load, respectively. The maximum magnetic flux penetrating the disk is defined by 

(^) max = f° B z 2n (wp/VK) dr : (13) 

where the upper limit in the integral corresponds to ^ out — > oo. According to assumption (iv), we have 
the two possibilities. 

(1) The parameter ^ out is finite with \p£, = and 6m = for ^/h < (^D)max- 

(2) The parameter ^ out is infinite with > and 6m > for ^/h > {^Djmax- 

It is easily to check that ( x I / D) max decreases with the increasing a* and n by using equations (8) 
and (13). Thus 6m occurs with infinite ^ out as a* and n are greater than some critical values. The 
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existence of the positive 9m is just the condition for CEBZMC, which can be expressed more clearly by 
the contours of 9 M in a* — n parameter space as shown in Figure 3. 

In Figure 3 the value of 9 m is labeled beside each contour, increasing along the direction from left- 
bottom to right-top, and the thick solid line labeled zero is the critical contour with 9 M = 0. Thus the 
parameter space is divided by the critical contour into two parts: (i) the right-top part for CEBZMC, 
and (ii) the left-bottom part for MC process only. The following conclusions can be obtained from 
Figure 3: 

(i) The state of CEBZMC depends on both the BH spin a* and the index n, and it occurs only if 
a* and n are great enough to fall in the right-top part of the parameter space; 

(ii) The greater is a*, the less is n for the given value of 9m- The greater are a* and n, the greater 
is 9m for the given value of the lower boundary angle 9^; 

(iii) The greater is 9l, the greater is the part for CEBZMC in the parameter space. 

The requirement of CEBZMC for a greater n implies that the magnetic field is more concentrated 
in the central region of the disk by equation (8), although we might have ^ out approach infinity in this 
case. 



Since the load disk in MC process is much better understood than the remote load in the BZ 
process, we can obtain a better picture of the transfer of energy and angular momentum between the 
BH and the disk. In this section we are going to give two equivalent descriptions on the transfer of 
energy and angular momentum in MC process. 



In WLM the transfer direction of energy and angular momentum in MC process is discussed in 
detail in terms of co-rotation radius r c , and it is defined as the radius on the disk, where the angular 
velocity of the disk, fl D , is equal to the angular velocity of the BH, Q H . The radius r c can be determined 
from the following equation: 



3. Transfer of energy and angular momentum 



3.1. Description by co-rotation radius 




-i -l 



1, 



(14) 



and is expressed by 



£c (a*) 



r c /r, 



ms 



x^(i-<?r 1/3 (i+9), 



(15) 



where q = a/1 — a* is a function of a*. The parameter £ c decreases monotonically with a* as shown by 
the dot-dashed line in Figure 4. By using equation (12) with assumption (iv) we have the curves of £ out 
versus a* for different values of the index n as shown also in Figure 4. 
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It is found that the MC region is divided into two parts by £ c : the inner MC region (henceforth 
IMCR) for 1 < £ < £ c and the outer MC region (henceforth OMCR) for £ c < £ < £ out . Therefore energy 
and angular momentum are always transferred by the closed magnetic field lines from the BH to the 
disk in OMCR with il D < Vl H , while the transfer direction reverses in IMCR with Q D > Q H . The 
correlation of the BH spin with the transfer direction is given as follows: 

(i) For 0.3594 < a* < 1 we have £ c within the inner edge, and the MC region is all in OMCR with 
the transfer direction from BH to disk. 

(ii) For a°"' < a* < 0.3594 we have 1 < £ c < ^ out , and the transfer is bi-directional, i.e. it is either 
from BH to OMCR or from IMCR to BH. 

(iii) For < a* < a°"' we have £ c > £ out , and the MC region is all in IMCR with the transfer 
direction from disk to BH. 



3.2. Description by flow of electromagnetic angular momentum and redshifted energy 

In MT the poloidal components of the flow of electromagnetic angular momentum and redshifted 
energy in the magnetosphere are expressed as follows: 

S P L = - (B ■ tu/At^ B p = (I/2nac) B p , (16) 

S P E =(S P E ) 1+ (S P E ) 2 , (17) 

where / in equation (16) is the current flowing downward through an m-loop, which is related to the 
toroidal magnetic field B by Ampere's law: 

B T = -m, 18 

azu z c 

where m is a toroidal Killing vector as given in MT. The two terms on RHS of equation (18) are 
expressed by 

(4) i = (ac/Air) (#» x B T ) , (s e ) ^ = uS L , (19) 

— * 

where E p in equation (19) is the poloidal electric field and satisfies 

E p = - (v F /c) x B p (20) 

In absolute-space formulation for the BZ process all the laws of physics are expressed in terms 
of physical quantities measured by "zero-angular-momentum observers" (ZAMOs) (Bardeen, Press & 
Teukolsky, 1972), and ZAMO angular velocity u is expressed by 



- 8- 



u = 2aMrjY?. (21) 
In equation (20) velocity v F of the magnetic field line relative to ZAMO is expressed as 

v F = a- 1 (n F -u)m (22) 
Incorporating equations (17) — (22), we have the expression for as follows: 

S F = Q F (I/2nac) B p (23) 

According to assumption (i) equations (16) and (23) are also applicable to MC process. Inspecting 
equations (16) and (23) we have the following results. 

(i) There is neither angular momentum flow nor energy flow at all unless poloidal currents are 
presented (7^0). 

(ii) Both angular momentum flow and energy flow along the magnetic field lines away from the BH 
to the disk if the current I flows downwards through an m-loop, i.e., I > 0, while the direction of the 
two will be reversed, if the current I flows upwards through an m-loop, i.e., I < 0. 

If the description by S P L and S P E is in accord with that by co-rotation radius r c , the sign of the 
current I should be dependent on whether the radius r is outside or inside r c . We shall discuss the sign 
of the current I by using the equivalent circuit in a unified model for the BZ and MC processes given 
in WXL as shown in Figure 5, which is adapted to the configuration of the magnetic field in Figure 1. 

In the i-loop of Figure 5 segments PS (characterized by the flux \&;) and QR (characterized by the 
flux + A\l/j) represent two adjacent magnetic surfaces, and segments PQ and RS represent the BH 
horizon and the load sandwiched by these two surfaces. AZm and AZu are the corresponding resistance 
of the BH horizon and the load, respectively. Asm = (A^i/2n)n H and Ae Li = — (A^i/2-n) Vlf are 
electromotive forces due to the rotation of the BH and the load, respectively. The minus sign in the 
expression of Ae Li arises from the direction of the flux. The current in each loop is expressed as 



Ae m + Aeu a" - (If 

h = Az Hi + Az Li = (A *- /2 ' r) (24) 

Suppose that an m-loop is embedded the magnetic surface QR. Inspecting Figures 1 with 5, we 
find that the current / flowing downwards through the m-loop is simply i.e., 



/ = (/!- h) + (J 2 - I 2 ) + ■ ■ ■ + (/,_! - +Ii = Ii (25) 

Suppose that the magnetic surface QR intersects with the disk at the co-rotation radius r c , we have 
f2f = il H at Ti = r c . Henceforth the magnetic surface QR is referred to as co-rotation magnetic surface 
(CRMS). Incorporating Figures 1 and 5 and equations (24) and (25), we conclude the following results: 
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(i) The mag netic surface inside CRMS penetrates OMCR with Of = Of < £l H and / > 0. Thus 
S p and S P E are in the same direction as B p , i.e., both angular momentum and energy are transferred 
from the BH to disk. 

(ii) The magnetic surface outside CRMS penetrates IMCR with Qf = Of > tt H and I < 0. Thus 

— * — * — * 

S p and S P E are opposite to B , i.e., both angular momentum and energy are transferred from the disk 
to the BH. 

Thus we have shown that the two descriptions for transfer of angular momentum and energy in 
MC process are equivalent. 



3.3. Equivalence of the two descriptions above and below the equatorial plane of the BH 

We can also prove that the two descriptions are equivalent both above and below the equatorial 
plane of the BH by considering the following points. 

(i) Direction of B p above and below the equatorial plane of the BH; 

(ii) Rotation of the BH with respect to the magnetic field lines frozen in the disk; 

(iii) Dragging effect of the rotating BH on the magnetic field lines; 

(iv) Equation (18) relating the direction of B T and the sign of the current / through an m-loop 
(downwards or upwards). 

We can determine the direction of B T by points (ii) and (iii), and determine the sign of / by 
point (iv). Combining the above points with the symmetry above and below the equatorial plane, we 
summarize the direction of S p and S P E as shown in Table 1. It is shown in Table 1 the direction of S p 
and S P E depends only on the place where the field line penetrates the disk (in OMCR or IMCR), and 
the two descriptions are equivalent both above and below the equatorial plane of the BH. 



3.4. A further discussion on flow of energy in MC process 

From equation (17) we find that S P E consists of two terms, {s p E ^j and {s E ^j , where the latter has 
the same direction as S E , and the former is related to the directions of E p and B T . We are going to 
give a discussion on the direction of {s p E ^j in detail. 

— * — * — * 

Direction of E p : From equation (20) we know that the direction of E p is normal to B p , depending 
on the cross product of v F and B p . While v F depends on the difference between Q F and ZAMO angular 
velocity uj by equation (22). Since the closed magnetic field lines are frozen in the disk, so Q F depends 
on the place where the field line penetrates the disk, i.e., 

n F = n D = _i . (26) 

M (£ 3/2 xL + a,) 1 } 
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From equation (21) we know that u depends on ZAMO coordinates (r, 9) in BH magnetosphere, 
and it satisfies VL H > u, where the equality holds as ZAMOs reach the horizon. Generally speaking, the 
farther are ZAMOs away from the horizon, the less is uo. Incorporating equations (20), (21), (22) and 

— * 

(26), we find that the direction of E p is a little complicated, depending on the position of ZAMOs at 
the field line and on the position of the foot of the field line on the disk. 

Direction of B T : The direction of B T has been shown in Table 1. Considering the dragging effect 
of the rotating BH on the magnetic field, the direction of B T depends on the difference between Q H 
and il F . The toroidal magnetic field B T is in the same direction as vector m for £l H < f2 F , and it is 
opposite to m for Vl H > il F . 



Based on the above discussion we can determine the direction of yS E j and those of the concerning 

quantities according to the three possibilities for the sequence in magnitude of Q H , u and Q F as listed 
in Table 2. It is found from Table 1 and Table 2 that ( S E ^ is in the same direction as S E in cases A, 

are m 



B, D and E, while it is opposite to S E in cases C and F. Considering that both S E and (^e) 
the same direction, we infer that (s e ^ is dominated by (#e) m cases C and F. 



2 



4. Profile of electric current in BH magnetosphere 



Based on MT and our model depicted in Figures 1 and 5 we can discuss the profile of the current 
on the BH horizon and that of the current density flowing onto the horizon from BH magnetosphere. 
As argued in MT the region outside the disk and the BH can be regarded as a force-free region, and 
the following equation is satisfied: 



pj+fj/c") xB = (27) 
According to Assumption (i) equation (27) is applicable to MC process, from which we infer that 

— » — * — *„ 

the poloidal component of the current density j p is parallel to B p by considering E 1 = in the stationary 
magnetosphere. We can derive the profile of the current on the horizon and the profile of the current 
density jmh in terms of the latitude of the horizon. 

The equivalent circuit in Figure 5 is applicable to both the BZ and MC processes, where the 
current in each loop circuit is exactly the current on the horizon. Substituting = Bh^^AI, 
AZh = IpAQjw and the concerning parameters in the Kerr metric into equation (24), we have the 
current on the horizon due to the BZ and MC processes as follows: 



Ihbz - B « M 2csc 2e_ {1 _ q) - /o 2csc^-(l-g)' (28) 
T - R M a * ^ ~ ^ - T a * (1 - S> 

Ihmc - B » M 2csc 2 6 _ {1 _ q) ~ J o 2 csc^-(l-g)' (29) 
where J = B H (GM/c) w 1.48 x lO 10 ^ (M/M ) A, and B 4 is the strength of the magnetic field on 
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the horizon in the unit of 10 4 gauss. The two parameters k = Vt F /Vt H and (3 = Vt F /Vt H = Vt D /VL H 
are the ratios of the angular velocity of the field lines to that of the BH in the BZ and MC processes, 
respectively. The value range of k is < k < 1, which cannot be determined precisely due to lack of 
knowledge of the remote load. Although people know very few about the remote load in the BZ process, 
the ratio k is assumed to be 0.5 for the optimal BZ power (MT). Compared with the BZ process the 
value of (5 can be determined precisely in terms of a* and £ by the following equation: 

(30) 

Incorporating the new mapping relation (11) and equations (29) and (30) we have the curves of 
Ihmc/Io versus 9 from 9 M to 9 L = 0.457T for the given values of the power law index n as shown in 
Figure 6. 

From Figure 6 we find the following results: 

(i) The current Ihmc varies with 9 on the horizon, where Ihmc/Io > indicates the current flowing 
from high to low latitude, and Ihmc/Io < indicates the current flowing from low to high latitude; 

(ii) From equation (29) we know that the sign of Ihmc/Io depends on the MC region with which 
its angular region is connected. As argued in Subsection 3.1, we have Ihmc/Io > for a% < a* < 1 
with the angular region corresponding to OMCR as shown in Figure 6a, and have Ihmc/Io < for 
< a* < with the angular region corresponding to IMCR as shown in Figure 6c. Both Ihmc/Io > 
and Ihmc/Io < occur for a% < a* < 0.3594, since both OMCR and IMCR exist in this case as shown 
in Figure 6b. 

From the above discussion, the current Ihbz/Iq exists only in CEBZMC with 9 M > 0. An inter- 
esting issue is whether the current is continuous at the boundary angle 9 M on the horizon. Based on 
the condition for CEBZMC we know that 9 M corresponds to ^ out — > oo, so we have f3 — > at 9 M by 
equation (30) except the Schwarzschild BH with a* = 0. Thus the difference between Ihmc/Io and 
Ihbz/Io is given by 

Um/Io = oq—^ — T\ v ( 31 ) 

2 esc 2 0m - (1 - q) 

It is found from equation (31) that Iqm > will hold for CEBZMC with 9m > 0, and the current 
Iqm flows from the magnetosphere onto the BH horizon, which could be regarded as a beam of electrons 
emitting from the horizon to the magnetosphere. Inspecting equations (28) and (29), we find that the 
current Iqm arises from the difference between the two parameters k and j3 at 9m, which correspond to 
the two different kinds of loads. 

Still by the conservation of current we can calculate the current Imh flowing from the magnetosphere 
onto the horizon in the range 9m < 9 < 9l, since Ihmc is not constant on the horizon. Inspecting 
Figure 5, we know that Imh arises from the difference between the current of the two adjacent loops. 
Incorporating the new mapping relation (11) and equations (29) and (30), we have the current density 
Jmh flowing from the magnetosphere onto the horizon as follows: 



(3 = Q D /Q H 



2 (1 + q) 
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JMH — 



1 dJ 



2-sm 2 0(l-<?) 
-1 



4(1-/3) cosg _ . ndp] 
2-sin 2 0(l- g ) blU 17 



(32) 



where j = #£tt = 0.108 x B 4 (M / ' M Q ) A ■ cm- 2 , and 



d£ = 3 (1 + q) X 3 ms e /2 sing 
^ a*(e 3/2 xL + «*) 2G ( a *^' n ) 
By using equations (32) and (33) we have the curve of JmhUo versus 6 as shown in Figure 7. 

Comparing Figure 6 with Figure 7, we find that the positive value of Jmh/ 'jo always accompanies 
with the increasing Irmc/Iq f° r the given values of a* and n, while the negative value of Jmh/ 'jo always 
occurs simultaneously with the decreasing Ihmc/Iq- These results are just required by the conservation 
of the current in BH magnetosphere. 



5. MC process and disk radiation 



Based on the three conservation laws of mass, energy and angular momentum, Li derived the 
radiation flux due to the MC effects as follows (Li02a): 



F MC = {E + - n D L+y 2 ^ (E + - n D L + ) Hrdr, (34) 

where H is the angular momentum flux transferred between the BH and the disk, and is assumed to be 
distributed from r = r ms to r = r out with a power law: 



tt J i T'ms *^ T *^ 1 out /qr-i 

t 0, r > r out ^ 

where A is regarded as a constant in Li02b. It is shown in Li02b that the magnetic coupling between a 
black hole and a disk can produce a very steep emissivity with index a = 4.3 ~ 5.0, where the emissivity 
index is defined as 



a = -dlnF/dlar. (36) 

This result is consistent with the recent XMM-Newton observation of the nearby bright Seyfert 1 
galaxy MCG-6-30-15, being regarded as one of the observational signatures of the magnetic coupling 
between a rotating BH and the surrounding disk. 

Unfortunately, the derivation of equation (35) was not given in Li02a and Li02b, which should be 
related to the torque exerted on the BH in MC process. In WXL the torque exerted on the BH is 
expressed by 
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. 9l (1 - B)sin 3 9d9 



where 



T = 5|M| x 3.26 x Wg ■ cm 1 ■ s~ 2 , (38) 

and M 8 is the BH mass in the unit of 1O 8 M . The flux H can be derived from T M c by the conservation 
of angular momentum as follows: 

aw* - - , (1 y> ( \~ T'T (39 > 

2 - (1 - q) sm 6 or 
where 89 jdr is related to the new mapping relation (11) by 

89/dr = (89/80 (dU8r) = _ G ( a *'^>™\ ( 40 ) 

r,„, >h\0 



Incorporating equations (39) and (40), we have 

H (a,; n)/H = { A r "' 1 (41) 



where H = 1.48 x 10 21 x B%M S g ■ s' 2 , and 



A (a*, £) - 2 jocsc^-n-o^i-^ (a*, 



^a(o*,£) 



27r X 2 ns [2csc 2 e-(l-g)] J 
V /(l+^Xm 4 s +2a2 X m 6 3 )(l-2Xm 2 s «- 1 +a2 Xm 4 s 5- 2 ) 



(42) 



Thus we derive the expression (41) for the function H, where m = —n, and the coefficient A (a*, £) 
is dependent on a* and £ rather than a constant given in equation (35). 

Incorporating equations (34), (36) and (41), we have the curves of the emissivity index a versus 
lg (r/M) for the different values of 9l and n in Figure 8. As shown in the shaded region of Figure 8, 
we find that the recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 
can be simulated by the new mapping relation (11) with three parameters: (i) the BH spin a*, (ii) the 
power law index n, and (iii) the lower boundary angle 9l- While the index produced by a standard 
accretion disc (SAD) is far below the shaded region as shown by the solid lines in Figure 8. More 
parameters suitable to the observation are listed in Table 3, from which we find that the condition for 
CEBZMC is well satisfied by these values of the parameters a* and n. According to the condition for 
CEBZMC derived in our model we expect that an jet might be produced nearby bright Seyfert 1 galaxy 
MCG-6-30-15 by the BZ process accompanying with MC process. 
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6. Summary 

In this paper the MC effects are discussed based on a stationary, axisymmetric magnetosphere 
anchored in the BH and its surrounding disk. Since we know very few about the magnetic field connecting 
BH with disk, the discussion is based on some simplified assumptions. It is known that electromagnetic 
field is related to electric charge density and current density in BH magnetosphere by a 'stream equation' 
containing the 'stream function' proposed in MT. Unfortunately it is a rather complicated task to solve 
analytically the 'stream equation' with the boundary conditions at horizon and disk surface. 

Recently the magnetic connection of a Kerr BH with disk and the resulting transportation of 
energy and angular momentum are discussed analytically based on a toy model in Li02c, where the 
poloidal magnetic field connecting the BH with the disk is naturally produced by a single toroidal 
current flowing around the BH in the equatorial plane. Although the existence of the single toroidal 
current needs further explanation, it is a step towards the goal of the origin of the magnetic connection. 
Although the configuration of the magnetic field in Li02c is different from that described in our model, 
the very steep emissivities in the disk are produced in both models. So we think that the very steep 
emissivities can be regarded as the main feature of the magnetic connection between the BH and the 
disk. 

The main results of our model are summarized as follows, (i) A consistent picture about the 
configuration of the magnetic field is depicted by Figure 1 with the equivalent circuit in Figure 5, and 
the transfer of energy and angular momentum along the closed field lines in MC process is described 
clearly and virtually in the two equivalent descriptions, (ii) The condition for CEBZMC is derived 
naturally based on some reasonable assumptions, which might be helpful in explaining the high energy 
radiation from BH-disk systems by combining the two mechanisms, (iii) The profile of the current on 
the horizon is calculated, and its continuity at the boundary angle is discussed based on the equivalent 
circuit, (iv) The MC effects on the disk radiation and the emissivity index are investigated, and the 
results turn out to be in accord with the observation. 

In our model the most important parameters of the system are the BH spin a* and the power law 
index n. The importance of BH spin a* rests in the two aspects: (i) The ratio of the rotating energy to 
the total energy of a Kerr BH is a function of a*, and it increases monotonically with the increasing a* 
(Wald 1984). (ii) The extracting powers in the BZ and MC processes are all proportional to al (WXL). 
The power law index n plays an important role in adjusting the profile of the magnetic field in the two 
aspects: (i) The more is the value of n, the more is the magnetic field concentrated in the central region 
of the disk; (ii) The more is the value of n, the less is the maximum magnetic flux penetrating the disk, 
and the more is the boundary angle 9 M , and the more is the BZ power due to the effects of the open 
field lines on the horizon. 

In this paper, we fail to discuss the MC effects on disk accretion rate by the following consideration: 
The variation of disk accretion rate due to the transfer of angular momentum is probably a dynamical 
process rather than a stationary one, which is much more complicated than the case we discuss in this 
paper. We shall discuss the MC effects on accretion rate in the future work. 
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Fig. 1. — The poloidal magnetic field connecting a rotating BH with remote astrophysical load and a 
surrounding disk 
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Fig. 2. — The curves of 6m (solid line) and ^ out (dashed line) versus a* with 6 L = 0.457T and different 
values of n: (a) n = 3.0, (b) n = 4.0, (c) n = 4.5. 
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Fig. 3. — The contours of constant values of 9 M in a* — n parameter space with (a) 9 L = 0.457T, (b) 
6 L = 0.40tt. 
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Fig. 4. — The curves of £ c (dot-dashed line) and ^ out versus a* for < a* < 1 with n — 1.1, 3.0 and 4.0 
in solid, dashed and dotted lines, respectively, (a) 9 L = 0.457T, (b) 9 L = 0.407T. 
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Ae Ln AZ Ln 



Fig. 5. — An equivalent circuit for a unified model for the BZ and MC processes 
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Fig. 6. — The curves of the current Ihmc/ % versus from 6*m to 6l = 0.457T with n = 3.0, 4.0 and 4.5 
in solid, dashed and dotted lines, respectively. 
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Fig. 7. — The curves of Jmh/jo versus 6 from Q M to ^ with n = 3.0,4.0 and 4.5 in solid, dashed and 
dotted lines, respectively, (a) a* = 0.998, (b) a* = 0.3, (c) a* = 0.1. 
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Fig. 8. — The emissivity index of the Seyfert 1 galaxy MCG-6-30-15 inferred from the observation of 
XMM-Newton is shown in the shaded region. The emissivity index of a standard accretion disk (SAD) is 
plotted in solid line. The emissivity index produced by MC process is plotted versus lg(r/M) in dotted, 
dashed and dot-dashed lines for n = 6.0,7.5 and 9.2 in (a), n = 4.3, 6.0 and 7.3 in (b), n = 6.8,8.0 
and 9.5 in (c), and n = 5.8,7.0 and 7.8 in (d), respectively. The following parameters are imposed: 
6 L = 0.40vr in (a) and (b), 6 L = 0.45tt in (c) and (d); a* = 0.995 in (a) and (c), a* = 0.998 in (b) and 
(d). 
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Table 1. Direction of S P L , S P E and corresponding quantities above and below EQP of the BH (EQP in 
Table 1 and Table 2 below is an abbreviation for equatorial plane) 



Position B p B T Current through m-loop S£ and Sg 



above 


OMCR 


From the BH 


Opposite to fa 


I > 0, downwards 


From the BH to the disk 


EQP 




to the disk 


















IMCR 




Same as fa 


/ < 0, upwards 


From the disk to the BH 



. . OMCR _ ... Same as fa I < 0, upwards From the BH to the disk 

below Mom the disk 

EQP to the BH 



IMCR Opposite to fa I > 0, downwards From the disk to the BH 



Table 2. Direction of (S E )i and corresponding quantities above and below EQP of the BH 





Case 




V F 




B T 






A 


n F > n H > w 


Same as fa 


Outward and normal to B p 


Same as fa 


Opposite to B p 


above 














EQP 


B 


n H >n F >ui 


Same as fa 


Outward and normal to B p 


Opposite to fa 


Same as B p 




C 


n H > w > n F 


Opposite to fa 


Inward and normal to B p 


Opposite to fa 


Opposite to B p 


below 
EQP 


D 
E 
F 


n F > n H > u> 
n H > n F > u> 
n H > oj > n F 


Same as fa 
Same as fa 
Opposite to fa 


Outward and normal to B p 
Outward and normal to B p 
Inward and normal to B p 


Opposite to fa 
Same as fa 
Same as fa 


Same as B p 
Opposite to B p 
Same as B p 



Table 3. The parameters a*, n, and 9 L adapting the emissivity index to the observations 



9 L a» = 0.990 a* = 0.992 a* = 0.994 a» = 0.996 a» = 0.998 



0.40tt 5.9- 11.1 5.6 - 10.4 4.9 - 9.6 3.8 -8.7 4.0- 7.3 
0.45tt 6.4- 11.3 6.1 - 10.6 5.8 - 9.9 5.3 -9.0 4.1 - 7.8 
0.48tt 6.5 - 11.4 6.3- 10.7 6.0- 10.0 5.6 -9.2 4.9-8.0 



